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Abstract 

Magnetoresistance (MR) of ion irradiated monolayer graphene samples with 
variable-range hopping (VRH) mechanism of conductivity was measured at 
temperatures down to T — 1.8 K in magnetic holds up to B = 8 T. It was ob¬ 
served that in perpendicular magnetic fields, hopping resistivity R decreases, 
which corresponds to negative MR (NMR), while parallel magnetic held re¬ 
sults in positive MR (PMR) at low temperatures. NMR is explained on the 
basis of the orbital model in which perpendicular magnetic held suppresses 
the destructive interference of many paths through the intermediate sites 
in the total probability of the long-distance tunneling in the VRH regime. 
At low helds, a quadratic dependence (\AR/R\ ~ B 2 ) of NMR is observed, 
while at B > B*, the quadratic dependence is replaced by the linear one. It 
was found that all NMR curves for different samples and different tempera¬ 
tures could be merged into common dependence when plotted as a function 
of B / B*. It is shown that B* ~ T 1//2 in agreement with predictions of the 
orbital model. The obtained values of B* allowed also to estimate the lo¬ 
calization radius £ of charge carriers for samples with different degree of 
disorder. PMR in parallel magnetic helds is explained by suppression of 
hopping transitions via double occupied states due to alignment of electron 
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1. Introduction 


Defects are very useful tool to modify and control the physical proper¬ 
ties of true two-dimensional (2d) material - monolayer graphene, therefore, 
disordered graphene attracts a lot of attention 0,3 3 8, i, 3 0. The evo¬ 
lution of optical (Raman spectra) and electronic (conductivity) properties 
of graphene with increasing disorder has been investigated in many papers 
(see, for example, [i, @, m Gj])- It has been shown in [11] that increase 
of disorder leads to gradual change of the mechanism of conductivity from 
metallic one in pristine films to the regime of weak localization, and then to 
the mechanism of variable-range hopping (VRH) of strongly localized charge 
carriers. In this work, we report the results of magnetoresistance (MR) mea¬ 
surements in samples of monolayer graphene strongly disordered with ion 
irradiation. MR in the VRH regime of conductivity has been observed ear¬ 
lier in different graphene-based structures: fluorinated graphene 0], graphene 
exposed to ozone [5j, graphene oxide 0. In Ref. [7], MR was measured in 
monolayer graphene flakes subjected to Ga + ion irradiation. In one highly 
disordered sample, the negative MR was observed which was attributed to 
the crystalline-boundary scattering. In Ref. [9j, disorder in graphene was 
introduced by C + ion irradiation with energy 35 keV. It was shown that at 
high dose of irradiation, the conductivity is described by the VRH mecha¬ 
nism, but no MR measurements were conducted. We are not aware about 
study of hopping MR in series of monolayer graphene samples with gradually 
increased degree of disorder. 


2. Experimental results and discussion 

The investigated samples belong to a series of samples disordered by 
different dose of ion irradiation. All samples were fabricated by means of 
electron-beam lithography (EBL) on the common large-scale (5x5 mm) 
monolayer graphene him and divided into 6 groups. Initial sample before 
EBL was marked as sample 0. The samples from the first group, marked 
as sample 1, were not irradiated, while 5 others groups were subjected to 
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different doses (from 5 x 10 13 up to 1 x 10 15 cm 2 ) of irradiation with carbon 
ions with energy 35 keV. In Ref. 10|], concentration of structural defects 


Np was determined for each group of samples using measurements of the 
Raman scattering. For samples 1-4, the values of N E in units of 10 12 cm -2 
were 1, 3, 6 and 12 correspondingly. In Ref. ll|, it was shown that the 


temperature dependence of resistivity R(T) in initial sample 0 has a metallic 
character, while in slightly irradiated sample 1, conductivity is characterized 
by the regime of ’’weak localization”. Measurements of MR in this sample 
showed a remarkable agreement with theoretical model (l2j |. In more irradi¬ 
ated samples 2, 3 and 4, dependence R{T) is described by the variable-range 
hopping (VRH) mechanism of conductivity typical for strongly localized car¬ 
riers [13]. As known, the form of R(T) in VRH depends on the structure of 
the density-of-localized states g(/i) in the vicinity of the Fermi level (FL) /i: 
when g(u) V 0, R(T ) is described by the ’’Mott VRH” which in 2d has the 
form of ’’T 1 / 3 - law”: 


R(T) = Ro exp (T m /T) 1 ^ 3 , T m = (1) 

Here i? 0 is a prefactor, Cm = 13.8 is the numerical coefficient, £ is the radius 
of localization. The Coulomb interaction between localized carriers leads to 
their redistribution in the vicinity of the FL. This results in a ’’soft” Coulomb 
gap around FL, which in 2d has a linear form g{e) — |e — g\(n/e 2 ) where k 
is the dielectric constant of the material. In this case, g{fi) = 0 and VRH is 
described by the Efros-Shklovskii (ES) VRH or ’’T 1 / 2 - law”: 

R(T) = Ro exp (Tes/T) 1 ' 2 , T es = C ES {e 2 /<) (2) 


Here C EE = 2.8 is the numerical coefficient. Measurements of MR in samples 
2, 3 and 4 was performed at temperatures down to 1.8 K in magnetic fields 
up to B — 8T in perpendicular and in-plane (parallel) B\\ geometry. It 
was observed that B± leads to negative magnetoresistance (NMR) while B\\ 
results in positive magnetoresistance (PMR) at low temperatures, Fig. 1. 
This anisotropy shows unambiguously that MR in perpendicular and parallel 
fields has different origin: NMR is determined by the orbital mechanisms, 
while PMR is determined by the spin polarization. In that order we will 
discuss the results of measurements. 


3. NMR in perpendicular magnetic fields 

Fig. 2 shows the MR curves AR(B)/R(0) = [R(B) — R(0)]/R(0) at 
different T for all three samples on a linear scale. One can see, that NMR at 
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Figure 1: A R/R as a function of B for sample 3 in parallel and perpendicular magnetic 
fields at different temperatures in K, shown near each curve. 


fixed T decreases with increase of disorder from sample 2 to 4. For samples 2 
and 3, NMR increases with decreasing T. For sample 4, A R/R first increases 
with decreasing T, but at T < 10 K, NMR rapidly decreases and the curves 
seek to change the sign. It could be due to the standard positive MR caused 
by the shrinkage of the wave functions in perpendicular magnetic fields 13 


Therefore we will discuss the NMR for sample 4 only down to 10 K. In Fig. 
3, NMR curves are plotted on the log-log scale. On this scale, the slope to 
the curve is equal to the power m in A R/R ~ B m . Quadratic dependence 
(m = 2) is observed at low fields up to some value B*. At B > B*, the 
quadratic dependence is replaced by linear one (m = 1) and then by sublinear 
dependencies. Some values of B* are shown in Fig. 3 by arrows. Very 
weak effect of NMR (about 1-2%) in VRH regime was earlier observed in 
three-dimensional (3d) conductivity in heavily doped and compensated Ge 
(for a review, see |14j). In 2d, a significant NMR in the VRH regime has 
been observed in perpendicular magnetic fields in different systems [15-20]. 
Anisotropy of this effect in perpendicular and parallel fields unambiguously 
indicates the orbital nature of NMR, because effects due to spin polarization 
has to be isotropic. Theoretically the effect of orbital NMR in the VRH 
regime of conductivity has been discussed in 2]J, 22, |23i, 24, [25|, (26j, |27j. The 
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main idea of the model suggested by Nguen, Spivak and Shklovskii 21] is 
based on the following consideration. 



-0.12 - 


Figure 2: NMR of samples 2 (a), 3 (b) and 4(c) at different temperatures (T, K) shown 
near the curves. The density of structural defects Nd (cm -2 ): in sample 2 - 3 x 10 12 , 3 - 
6 x 10 12 , 4 - 1.2 x 10 13 . 


In VRH, only part of localized states with energy levels within so-called 
’’optimal band” around FL e(T) are involved in the hopping process. In 
’’Mott VRH”, e(T) decreases with decrease of temperature [131 ]: 

e(T) = T 2 / 3 [g^)e}~ 1/3 (3) 

Correspondingly the hopping distance increases, which gives ~ T -1 / 3 : 

r h ~ « £(T m /T) 1/3 . (4) 

In ”ES VRH”, rh ~ T -1 / 2 . Therefore, at low T, rh becomes much larger 
than the mean distance between localized centers, and the probability of the 
long-distance hop is determined by the interference of many paths of the 
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Figure 3: NMR as a function of B plotted on the log-log scale for samples 2 (a), 3 (b) and 
4 (c). Short dashed lines correspond to m = 2, dashed lines to m = 1. B* indicate the 
end of quadratic dependence. 

tunneling through the intermediate sites which include scattering process 
(Fig. 4). All these scattered waves together with non-scattered direct wave 
contribute additively to the amplitude of the wave function \& 12 which reflects 
the probability for a charge carrier localized on site 1 to appear on site 2. 
There is no backscattering, and scattered waves decay exponentially with 
increasing distance as exp(2r/£), therefore only the shortest paths contribute 
to "012 • All these paths are concentrated in a cigar-shaped domain of the 
length Th, the width D fa (r/^) 1 / 2 and the area A ta r^ 2 £ 1//2 , where a < 1 is 
a numerical coefficient. 

As a result of averaging over different configurations, the contribution of 
the scattered sites to the total hopping probability vanishes due to destructive 
interference. The perpendicular magnetic held suppresses the interference 
which leads to the increase of the hopping probability and, therefore, to 
NMR. Taking into account that the hopping distance r^ is fixed, the increase 
of the hopping probability could be considered as a small increase of the 
localization radius £. Both ’’Mott” and ”ES” parameters Tm and Tes depend 
on the value of £, but differently: Tm £ 2 , Eq. (1), while T E s ~ Eq. 
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Figure 4: Schematics of the cigar-shaped region with localized states contributing to the 
probability of an electron tunneling from center 1 to center 2. 


(2). This gives rise to a conclusion that NMR has to be stronger for the case 
of ” Mott VRH”, then that for ” ES VRH”. This suggestion was observed in 
experiment. In Fig. 5(a), the temperature dependences of the modulus of 
NMR at fixed magnetic field B — 4 T for samples 2, 3 and 4 are shown. 
One can see that A R/R increases with decreasing T linearly with T^ 1 / 3 and 
deviates from the straight line to smaller values at approximately the same 
temperatures where initial curves R(T) at B = 0 deviate from T 1 /3 - law to 
the stronger T 1//2 - law, Fig. 5(b). 

In accordance with theoretical considerations 21, 22[ 23[ 24[ 25[ 26], NMR 
as a function B has to be linear at moderate fields and quadratic at very low 
fields. In the ” orbital” model, it is natural to normalize the magnetic field 
using the ratio r) = $b/^>o where <f>o = h/ 2e ~ 2.07 x 10~ 15 W is the magnetic 
flux quantum and = B ■ A is the magnetic flux through the cigar-shape 


area. We suggest that B* corresponds to rj = 1, i.e. B* = <&o/A 
Taking into account Eq. (4), one get (a ~ 1): 


A 3//2 (( -1 / 2 


B* = $ 0 C 2 (Tm/T)~ 1/2 = A T 1/2 , 


A = $oC 2 Tr 1/2 


■ M 


(5) 


In Fig. 6, the values of B* for all samples are plotted as a function of T 1 / 2 . 
One can see that, indeed, B* ~ T 1//2 . Coefficient A is equal to 0.1, 0.24 
and 0.58 [T K” 1 / 2 ] for samples 2, 3 and 4. Knowledge of A allows us to 
estimate with accuracy of a the values of localization radius £. For samples 
2, 3 and 4, with T M = 68, 308 and 5960 K jllj, this gives £ = 50, 22 and 7 
nm correspondingly, which is quite reasonable: £ decreases with increase of 
disorder (from sample 2 to sample 4). 

We also use the values of B* in an attempt to merge the NMR data for 
all samples and all temperatures below 25 K. In Fig. 7, NMR curves from 
Fig. 2 are plotted as a function of dimensionless parameter B/B*. One can 
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Figure 5: (a) modulus of NMR at fixed magnetic field B = 4 T for samples 2, 3, 4 plotted 
as a function of T -1 / 3 . (b) - temperature dependence R(T) for these samples at zero 
magnetic field. 


see that all curves are merged in a universal dependence. No saturation of 
NMR is observed up to rj > 60. 


4. Positive MR in parallel magnetic fields 

In parallel magnetic fields, PMR is observed at low temperatures, Fig. 1. 
Very small NMR at high temperatures could be explained as the traces of 
NMR due to the possible folds on the surface of monolayer graphene hint, 
where the parallel magnetic held has a perpendicular influence component. 
PMR appears only at low temperatures, at T < 5 K and increases with 
decreasing T. PMR is proportional to B y at low helds, and becomes linear 
with increasing B (Fig. 8). 

PMR in parallel magnetic helds has been observed earlier in 2d VRH 
regime in different systems: in Af,.In 1 _ a ,Sb/InSb quantum well 
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GaAs/Al x Gai_a,As heterostructure 27]. Because the parallel magnetic held 


m a 


couples only to the electron spin, it means that the spin state of localized 
electrons influences the hopping conductivity despite the fact that it is not 
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Figure 6: The values of B* as a function of T 1 / 2 . The sample numbers are shown near 
the straight lines. 


included explicitly in the expressions for VRH. The explanatory models of 
this effect are based on two ideas. The authors [28] considered the case when 
intermediate scattering centers (see Fig. 4) should be occupied to produce 
a negative scattering amplitude. Thus, the interference is depended on the 
mutual spin orientation of the hopping electron and electron localized on 
scattering center. In magnetic fields, all localized spins are aligned which in¬ 
creases the destructive interference and results in increase of the resistance. 
Another mechanism was suggested by Kurobe and Kamimura j29j and stud¬ 


ied in 30]. In this model, it is recognized that a certain fraction of the states 


can accommodate two electrons. Double occupancy is possible if the on-site 
Coulomb repulsion U between the electrons is smaller than the width of the 
energy distribution function of the localized states. It was already mentioned 
that in VRH, only localized states with energy level within the narrow op¬ 
timal band of width e(T) around /r are involved in the hopping process, Eq. 
(3). However, for some states, which cannot participate in VRH at given 
temperature because the energy of the Erst electron e^ 1 ) is well below /r, the 
energy of the second electron = e^ 1 ) + JJ may be located just within the 
optimal band, Fig. 9. This allows those states to participate in the VRH at 
zero magnetic field. In strong field limit, all spins are polarized and, there¬ 
fore, transitions through the double occupied states are suppressed which 
results in increase of resistance. 

Fig. 1 shows that PMR is observed in our samples only at low tem¬ 
peratures. We think that this fact supports the mechanism [29] based on 
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Figure 7: The NMR data for different samples and different temperatures plotted as a 
function of dimensionless magnetic field B±/B*. The arrows and numbers show the end 
of each curve and indicate the sample (2, 3, 4) and T. In insert, the NMR data are shown 
for small values B±/B* < 1. 


participation of the double occupied states, because the first mechanism 
[28[ has no limitation for observation at all temperatures. In the second 
mechanism, however, contribution of the double occupied states in VRH is 
important only when the width of the optimal band e(T) which decreases 
with decreasing T, becomes less than U, otherwise in the case of an oppo¬ 
site inequality, U e(T), the localized states will either participate or not 
participate in VRH independently of the existence of the double occupied 
states. At moderate magnetic fields, theory [30i] predicts the linear depen¬ 
dence A R/R ~ (g L g B B)/T, where is the Lande-factor and fi B is the 
Bohr magneton, while at weak fields one expect the quadratic dependence 
A R/R ~ B 2 . This agrees with experiment (Fig. 8). Theory predicts also 
saturation PMR at strong fields when all electron spins are polarized. In our 
samples, no tendency to saturation was observed in magnetic fields up to 8 
T. 


An interesting issue not discussed above is a possible role of spin-orbit 
coupling for magnetoresistance in disordered graphene in the VRH regime. 
The theory of anisotropic magnetoresistance of spin-orbit coupled carriers 
scattered from polarized magnetic impurities, in which spin-orbit coupling 
enters via specific spin-textures on the carrier Fermi surfaces and ferromag¬ 
netism via elastic scattering of carriers from polarized magnetic impurities, 


was proposed in Ref. 31]- Later, nonlinear anomalous Hall effect and nega¬ 


tive magnetoresistance in a system with random Rashba held was also studied 
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Figure 8: PMR of sample 3 plotted as a function of B\\. In insert, PMR data for low fields 
are plotted on a quadratic scale. 


theoretically [32|. There it was shown that electron scattering from a fluctu¬ 
ating Rashba held in a two-dimensional nonmagnetic electron system leads 
to a negative magnetoresistance arising solely due to spin-dependent effects. 

However, knowing relative weakness of spin-orbit interaction in graphene, 
we think that it can only slightly modify the interpretation of the presented 
results. 



Figure 9: Schematic representation of the possible hopping transition via the double oc¬ 
cupied states where the double occupied is an initial ei (a), final ej (b) or both initial and 
final states (c). /r represents the position of the Fermi level, dashed lines show the width 
e(T) of the optimal band at given temperature, Eq. (3). 
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5. Conclusions 

In conclusion, MR of strongly disordered monolayer graphene films with 
VRH mechanism of conductivity was measured in perpendicular 13_ and 
parallel B\\ magnetic fields. It is shown that in perpendicular B±, resistance 
decreases (negative MR), while B\\ leads to positive MR at low temperatures. 
The NMR effect is explained in the framework of the orbital mechanism 
based on the interference of many paths through the intermediate sites in 
the probability of the long-distant tunneling in the VRH regime. Magnetic 
field B* at which the quadratic dependence A R/R ~ B 2 is replaced by the 
linear one is determined as the merging parameter. It is shown that all 
MNR curves for different samples and different temperatures are merged in 
a universal dependence plotted as a function of dimensionless field B/B*. It 
was shown that B* ~ T 1//2 in accordance with the ’’orbital” model, the slope 
of this dependence allows us to estimate the localization radius £ for samples 
with different dose of irradiation. Expectedly, £ decreases with increase of 
disorder. The PMR effect in parallel fields could be explained by suppression 
in a strong magnetic field the hopping transitions via double occupied states 
due to electron spin polarization. 
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